Viscosity anomaly in core-softened liquids 
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The present article presents a molecular dynamics study of several anomalies of core-softened 
systems. It is well known that many core-softened liquids demonstrate diffusion anomaly. Usual 
intuition relates the diffusion coefficient to shear viscosity via Stockes-Einstein relation. However, 
it can break down at low temperature. In this respect it is important to see if viscosity also 
demonstrates anomalous behavior. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60. Kw 



00 



INTRODUCTION 



SYSTEMS AND METHODS 



o 



O 

o 



> 
o 

rn 
O 
m 



X 



It is well known that some liquids (for example, wa- 
ter, silica, silicon, carbon, and phosphorus) show an 
anomalous behavior [ll-Q:their phase diagrams have re- 
gions where a thermal expansion coefficient is negative 
(density anomaly), self-diffusivity increases upon com- 
pression (diffusion anomaly), and the structural order of 
the system decreases with increasing pressure (structural 
anomaly) A number of studies demonstrates water- 

like anomalies in fluids that interact through spherically 
symmetric potentials (see, for example, [3, (l^, |3| and 
references therein) . Many of these studies report the ap- 
pearance of diffusion anomaly in different systems. How- 
ever, the diffusion coefficient is closely related to shear 
viscosity of liquid therefore one can expect that the shear 
viscosity also demonstrates some kind of anomalous be- 
havior. 

Although many studies of core-softened systems report 
the diffusivity calculations there is a lack of studies which 
calculate shear viscosity. This can be related to the fact 
that viscosity is much harder to compute in simulation. 
So the usual intuition is applied: the viscosity can be ex- 
tracted from the diffusion coefficient by Stockes-Einstein 
(SE) relation 8]. However, it was recently found that 
SE relation can be violated at low temperatures [l3] ■ 
This case the usual intuition can fail to predict the vis- 
cosity behavior correctly. In this respect it is important 
to monitor both the diffusion coefficient and shear vis- 
cosity of corc-softcncd liquids at low temperatures to see 
theirs behavior in the regions of anomalous behavior. 

The goal of the present article is to investigate the be- 
havior of shear viscosity of core-softened fluids at low 
temperatures, to see if theirs shear viscosity demon- 
strates anomalous behavior and if so to find the relation 
between the viscosity anomaly region and the regions of 
other anomalies. 



Two systems are studied in the present work. The first 
one is a core-softened system introduced by de Oliveira 
et al TT] . This system is described by the Lennard- Jones 
potential with Gaussian well (LJG): 
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with a — 5.0, ro/cr — 0.7 and c = 1.0. The diffusivity of 
this system was studied in several papers [lll - [l5t . Note 
that the parameters of the potential are chosen in such a 
way that the effect of attraction becomes negligibly small 
and one can consider this system as a purely repulsive 
core-softened one. 

The second system studied in this work is Soft Repul- 
sive Shoulder System (SRSS) introduced in the work [1^ . 
The potential of this system has form: 
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where a is "hard"-core diameter, cti = 1.35 is soft-core 
diameter and fcg = 10.0. In Ref. JTJ it was shown that 
this system demonstrates anomalous behavior. Our later 
publications gave detailed study of diffusion, density and 
structural anomalies in this system [Ts! ]. 

It is well known that there is a close link between dif- 
fusion coefhcient and shear viscosity of liquids. Viscosity 
is a quantity which is usually measured in experiments. 
However, due to the technical problems the simulation 
of shear viscosity represented in the literature is very 
poor. One of the goals of this article is to study the 
behavior of viscosity of the two model liquids described 
above. Taking into account that the diffusion coefficient 
demonstrates anomalous behavior for these systems we 
are interesting to see if the viscosity also demonstrates 
some anomalies. 
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In order to study the transport coefRcients of the sys- 
tems we used Molecular Dynamics method. In both cases 
a system of N — 1000 particles was simulated. The equa- 
tions of motion were integrated by velocity Verlet algo- 
rithm. In case of LJG system the time step was set to 
dt = 0.001, the equilibration period was 3.5 • 10^ steps 
and the production period 1 • 10^ steps. In the case of 
SRSS the equilibration the time step was dt — 0.0005, 
the equilibration period was 3.5 • 10^ and the production 
run was 0.5-10^ steps. The cut-off radius was set 3.5 for 
LJG system and 2.2 for SRSS. Velocity rescaling was ap- 
plied during equilibration, the production corresponded 
to NVE ensemble. Shear viscosity is difficult to measure 
in simulation because of large fluctuations of shear stress 
function. In order to improve the precision of the data 
we increased the equilibration time in anomalous region 
up to 7.5 • 10^ and the production time up to 1.5 • 10^ for 
some simulations. 

In order to get good statistics on the transport proper- 
ties of the systems many data points were simulated. In 
the case of LJG system the data points were chosen in the 
density interval from p = 0.05 till p = 0.3 with step 5p = 
0.01 along several isotherms. The following isotherms 
were considered: T = 0.15; 0.2; 0.25; 0.3; 0.4; 0.5; 1.0. In 
order to see the anomalous region better we also simu- 
lated the isotherms T — 0.17 and 0.23 for the densities 
from p = 0.08 up to 0.18 with step 0.01. 

In case of SRSS we used the densities from p = 0.3 
up to p = 0.8 with step 5p = 0.05 and temperatures 
T = 0.2; 0.25; 0.3; 0.35; 0.4; 0.5; 0.7 and 1.0. 

The diffusion coefficients were computed via Einstein 
relation and shear viscosity by integration of shear stress 
autocorrelation function. 



RESULTS AND DISCUSSION 
Lennard-Jones - Gauss system 

As it was mentioned in the introduction the visc osity 
anomaly in LJG system was already reported in Ref. [l9j . 
Here we make a more detailed simulation study of this 
anomaly. Our goal is to see the location of the anomaly 
VD. p — T plane and its relation with other anomalies, 
such as diffusion anomaly, density anomaly and struc- 
tural anomaly. 

Figs. [T] shows the viscosity along several isotherms for 
LJG system. One can see that the anomaly is very pro- 
nounced for the temperature T = 0.15, but it rapidly 
disappears with increasing temperature. At T = 0.3 the 
anomaly is already of the order of numerical accuracy and 
we estimate this temperature as the temperature where 
viscosity anomaly disappears. 

The location of diffusion, density and structural 
anomalies of LJG system in p — T plane have already 
been reported in literature |12]. Figs. [5] shows the re- 
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FIG. 1: (Color online). Viscosity anomaly in JLG system. 



gions of these anomalies and the viscosity anomaly. In- 
terestingly, as it was proposed in Rcf. 22] the anomalous 
regions are enveloped in each other. However, the vis- 
cosity anomaly violates this rule: it has partial overlap 
with density anomaly, but no one of them is inside of one 
another. It was also shown in the literature that from 
the thermodynamic arguing it follows that the density 
anomaly region is always inside the structural anomaly 
one, while the diffusion anomaly can have any location 
with respect to the other anomalies [2^ . The viscos- 
ity anomaly is another example of anomalies of dynamic 
rather then thermodynamic properties. Therefore, one 
can expect that the viscosity anomaly can also have any 
possible location with respect to the density and struc- 
ture anomalies. 

In our previous work we showed that the anomalies 
can be visible along some paths in thermodynamic space 
while along others they can be invisible 13l4l5|. For ex- 
ample, diffusion anomaly is seen along isotherms but not 
isochors. An important consequence of this difference is 
that Rosenfeld excess entropy scaling for diffusion coeffi- 
cient 2J, |25| is fulfilled along isochors but breaks down 
along isotherms [l^. This makes important to see the 
viscosity behavior along different trajectories. 

Figs. |3] shows shear viscosity of LJG system along sev- 
eral isochors. One can see that viscosity demonstrates 
a minimum. Viscosity minimum along isobars was ob- 
served experimentally for water [2^. The authors called 
this minimum as "viscosity anomaly". However, in our 
previous work we showed that viscosity minimum along 
isochors appears naturally because of the interplay of 
potential-potential and kinetic-kinetic correlations even 
in simple liquids [2l| . The same results can be obtained 
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FIG. 2: (Color online). Location of anomalous regions in p—T 
plane for LJG system. 



FIG. 4: (Color online). Rosenfeld relation for shear viscosity 
of LJG system along isochors. 
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Fig. [5] (b) shows the Rosenfeld scaling for high tempera- 
ture (T = 1.0). Points correspond to the data from simu- 
lations while the straight line is the best fit line. One can 
see that the simulation points demonstrate some kind of 
oscillations around the best fit line. We relate these os- 
cillations to the numerical inaccuracies in the viscosity 
computations and we believe that Rosenfeld scaling of 
viscosity does work for high enough temperatures. 

One can conclude that as in the case of diffusion co- 
efficient Rosenfeld relation for shear viscosity in systems 
with thermodynamic anomalies holds true along isochors 
but breaks down along low temperature isotherms. This 
confirms the idea of different behavior of a system along 
different trajectories in p — T — P space which is discussed 



in details in Ref . [ij, ll5j . 



Soft Repulsive Shoulder system 



FIG. 3; (Color online). Shear viscosity of LJG system along 
several isochors. 



for isobars (not shown in Ref. [2l|). 

Figs, m shows the Rosenfeld relation for shear viscosity 
of LJG system along isochors. One can see that the linear 

— 2/3 

relation between ln( ''^iy2 ); which is predicted by Rosen- 
feld relation, holds true except the low Sex region. How- 
ever, if we consider the Rosenfeld relation along isotherms 
we see that it breaks down at low temperature (Figs. [5] 
(a)). Here we observe a self crossing loop like the one 
observed for diffusion in our previous works [isl - fisi [l8| . 



The second system considered in this work is Soft Re- 
pulsive Shoulder System (Eq. (2)). Phase diagram and 
anomalous behavior of this system were studied in de- 
tails in our previous works [16l4l8| . However, the shear 
viscosity of this system is measured for the first time. 

Figs. [B] (a) and (b) show the shear viscosity of SRSS 
system along several isotherms. One can see that at the 
lowest temperature T = 0.15 a tiny loop develops at the 
densities p = 0.40 — 0.45. However, the size of this loop 
in inside the error bar, so we can not consider it as a real 
anomaly. We believe that the anomaly appears at lower 
temperatures. Though the shear stress autocorrelation 
function decays very slowly and viscosity calculations be- 
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FIG. 5: (Color online). Rosenfeld relation for shear viscosity 
of LJG system along isotherms at low and high temperature. 



come very difficult. 

Figs. [7] shows the shear viscosity plotted along a set of 
isochors. One can see that for all presented densities the 
viscosity curves monotonically decrease with increasing 
temperatures. Basing on the arguments of our previous 
work |2l| we expect that shear viscosity passes a mini- 
mum at higher temperatures. 

Figs. |8] (a) and (b) show the Rosenfeld scaling plots 
for SRSS along isotherms and isochors. One can see that 
the scaling relations break down in case of isotherms. 
The reason for this breakdown is that while viscosity 
is monotonous function of density the excess entropy 
demonstrates anomalous behavior [l5| . As a result 
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FIG. 6: (Color online). Viscosity of SRSS system along a 
set of isotherms at (a) low and (b) intermediate and high 
temperatures. 



viscosity as function of Sex demonstrates nonlinear be- 
havior. 

At the same time Rosenfeld relation holds true along 
isochors (Figs. [8] (b)). All isochors can be divided in 
low(p < 0.45) and high density (p > 0.45) groups. The 
curves belonging to the same groups have similar slopes 
while the slope of the curves from different groups is es- 
sentially different. The reason for this change is that 
at low densities the system can be essentially approxi- 
mated by the system with effective diameter a while at 
high densities with diameter d. This change of particle 
size alters also the kinetic and thermodynamic proper- 
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FIG. 7: (Color online), 
of isochors. 



Viscosity of SRS system along a set 



ties of the system which we observe as the slope change 
in Figs.[S](b). 



Stokes-Einstein Relation 



In our previous work IJ, [15| we showed that the dis- 
crepancy in the diffusion and structural anomalies re- 
gions leads to the Rosenfeld relation breakdown along 
isotherms. As it was shown in the previous section the 
regions of diffusion and viscosity anomalies are also dif- 
ferent. In this respect it becomes important to see if 
Stockes-Einstein relation still holds true in the anoma- 
lous regions. 

The Stockes-Einstein relation can be written in the 
following form: 



CSE 



■nDrjd 



const, 



(3) 



where d is the character particle size. The coefficient 
cs E should be approximately constant and belong to the 
interval 2 < cg^; < 3. The limiting values cse = '2 
and csE — 3 correspond to the stick and slip boundary 
conditions. 

Figs, ini (a) and (b) show the Stockes-Einstein coeffi- 
cient for LJG and SRSS systems along isochors. One can 
see that in both cases cse is not constant. At the same 
time the numerical values of cse in both cases can be far 
from the interval [2, 3]. In our previous work we studied 
Stockes-Einstein relation for Soft Spheres system [2^ . It 
was found that Stockes-Einstein relation can be fulfilled 
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FIG. 8: (Color online). Rosenfeld scaling plots for SRS along 
(a) isotherms and (b) isochors. 



there if one takes the effective particle diameter as Barker 
perturbation theory one d - y-i/"^ where n is the soft- 
ness coefficient. However, in the case of core-softened sys- 
tems the situation is more complicated. These systems 
have two character length scales and the applicability of 
perturbation theory to such systems is questionable. 

From figs.|n](a) and (b) one can see that the qualitative 
behavior of cse is defined by the viscosity behavior. For 
example, the maximum of cse appears at the tempera- 
tures of minimum of viscosity. At the same time viscosity 
of SRSS is monotonous and we observe that cse is also 
monotonous in this case. However, one needs more de- 
tailed investigations to give any conclusive statements on 
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We find that the viscosity anomaly does exist, however, 
the region of (p — T) parameters where viscosity demon- 
strate anomalous behavior is different from the diffusion 
anomaly region. We place the regions of different anoma- 
lies in the same plot to see the relations between them. 
Finally we check the Stockes-Einstein relation for the liq- 
uids under investigation. 
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FIG. 9: (Color online), cse for (a) LJG system and (b) SRS 
along several isochores. 

Stockes-Einstein relation for core-softened systems. 

CONCLUSIONS 

It is well known from the literature that many core- 
softened liquids demonstrate some kind of anomalies. 
One of the typical anomalies in the core-softened systems 
is diffusion anomaly. It is also widely believed that diffu- 
sion is strongly connected to shear viscosity by Stockes- 
Einstein relation. In this respect it is interesting to see 
if the same systems demonstrate viscosity anomaly as 
well. In the present work we investigate this question. 
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